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THE PECULIAR X-RAY TRANSIENT ICR 16358 - 4726 
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ABSTRACT 

The new transient IGR 16358-4726 was discovered on 2003 March 19 with INTEGRAL. We detected 
the source serendipitously during our 2003 March 24 observation of SGR 1627 — 41 with the Chandra 
X-ray observatory at the 1.7 x 10 -10 ergs s _1 cm~ 2 flux level (2 — 10 keV) with a very high absorption 
column (Ah = 3.3 x 10 23 cm~ 2 ) and a hard power law spectrum of index 0.5(1). We discovered a 
very strong flux modulation with a period of 5880(50) s and peak-to-peak pulse fraction of 70(6) % 
(2 — 10 keV), clearly visible in the x-ray data. The nature of IGR 16358 — 4726 remains unresolved. 
The only neutron star systems known with similar spin periods are low luminosity persistent wind-fed 
pulsars; if this is a spin period, this transient is a new kind of object. If this is an orbital period, then 
the system could be a compact Low Mass X-ray Binary (LMXB). 
Subject headings: pulsars: individual (IGR 16358 — 4726) 



1. INTRODUCTION 

Ten new hard and relatively faint X-ray transients have 
been detected in the last few months during galactic 
plane scans by the IBIS/ISGRI detector (Ubertini et 
al. 2003; Lebrun et al. 2003) on the INTEGRAL Ob- 
servatory (Winkler et al. 2003). Their unusual spec- 
tral hardness has led to suggestions that these sources 
belong to a group of highly absorbed galactic binary 
systems, discovered due to the unique IBIS sensitivity 
in higher energy bands than earlier monitoring instru- 
ments. Two of these sources, IGR J16318 — 4848 and 
IGR J16320 — 4751, were subsequently observed with the 
AMM-Newton satellite. The EPIC-PN detector spec- 
trum of the former source revealed a very strong Fe 
K Q emission line and significant and Ni K a lines 
that were variable in time (Matt & Guainazzi 2003), 
whereas the latter exhibited no line features (Rodriguez 
et al. 2003). The spectral continuum of four out of the 
five sources (measured with the Rossi X-ray Timing Ex- 
plorer (RXTE) or AMM-Newton) was best fit with power 
law functions with photon indices, 7 ~ 1.6. Finally, 
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all sources showed variations ranging in timescales of 
~ 100 s to a few 1000 s, but none showed evidence for 
coherent modulations. 

IGR 16358 - 4726 was detected (Revnivtsev et 
al. 2003a) with INTEGRAL on 2003 March 19 during a 
standard (~ 2200 s) Galactic plane scan as a 50 mCrab 
source (15 — 40 keV). Preliminary analysis of additional 
INTEGRAL scans on 2003 March 22 and 26 indicated 
a decline in the source flux. Analysis of archival ob- 
servations with the Advanced Satellite for Cosmology & 
Astrophysics (ASCA) and BeppoSAX in 1999, revealed a 
source consistent with the location of IGR 16358 — 4726. 
BeppoSA X observations in 1998 & 2000, however, did not 
detect the source. 

We serendipitously observed IGR 16358 — 4726 dur- 
ing a Chandra scheduled observation of SGR 1627 — 41 
(Kouveliotou et al. 2003a), and then again with Direc- 
tor's Discretionary Time (DDT) one month later. In ad- 
dition we triggered Target of Opportunity observations 
with RXTE and with the Australia Telescope Compact 
Array (ATCA); the analysis of the RXTE data is re- 
ported in detail by Revnivtsev (2003c). ATCA observed 
IGR 16358 — 4726 over a 10.5 hr period beginning at 
UT 11:53 on 2003 April 13. The total on-source inte- 
gration time was 5 hr. No source was detected at the 
position of the X-ray source at an observing frequency 
of 4.8 GHz, with a 5cr upper imit of 0.9 mJy. We report 
in Sections 2, 3 and 4 the results of our analysis of the 
archival and current X-ray observations and discuss in 
Section 5 the implications of these results for the nature 
of IGR 16358 - 4726. 

2. ARCHIVAL X-RAY OBSERVATIONS OF IGR 16358 - 4726 

The field containing IGR 16358 — 4726 was observed 
with ASCA on 1999 February 26-28 during a scheduled 
observation of SGR 1627 — 41. The transient was only 
in the (wider) field of view of the two ASCA/Gas Imag- 
ing Spectrometers (GIS). We combined the GIS data and 
grouped the resulting spectrum in energy bins that con- 
tained at least 25 events each. Using XSPEC (vll.2; 
Arnaud 1996) we fit an absorbed single power law (PL) 
model, which gave an acceptable fit (x 2 jv = 56.1/58). 
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The best fit column density and photon power law in- 
dex are, 7V H = x 10 22 cm" 2 and 7 = 1.3±g;l, re- 
spectively. We derive an unabsorbed flux of 2.4 x 10~ 12 
ergs cm~ 2 s (2—10 keV), consistent with the analysis of 
Revnivtsev et al. (2003b). 



Table 1 

ASCA & BeppoSAX Observations of IGR 16358 - 4726. 
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06 Aug 1998" 45/ 21 <6 <1~3 
16 Sep 1998 6 30/ 12 < 9 < 21 
08 Aug 1999 34/ 80 10(2) 22(4) 
26 Sep 2000^ 61/ 21 <3 < 7 

a Flux is in units of 10~ 13 ergs s" 1 cm" 2 , b 90% upper limit 

BeppoSAX observed the SGR 1627-41 field four times 
(Table 1). Using XSELECT (v2.1) we extracted spectra 
from circular regions centered on IGR 16358 — 4726 with 
radii of 4' and 6' for the LECS and MECS, respectively. 
We use the standard response files provided by the Bep- 
poSAX Data Archive at HE AS ARC. The source was not 
detected during the observations in 1998 or 2000 and we 
derived a 90% upper limit for the unabsorbed flux (Table 
1). We detect the source during the 1999 observations, 
albeit very faint. We have used the spectral parameters 
derived with the ASCA data to fit the BeppoSAX source 
counts and find a flux value that is slightly higher than 
the one reported by Revnivtsev et al. (2003b). In fact, 
the true flux may be even higher, when one takes into 
account the partial source blockage due to the support 
structure for the MECS detector windows. 

3. Chandra OBSERVATIONS OF IGR 16358 - 4726 

We observed the transient for 25.7 ks with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) on Chan- 
dra on 2003 March 25 (Obsl). We used the CIAO 
v2.3 and CALDB v2.21 software for all data analysis 
and processing tasks. Details of the observation and 
data reduction techniques are presented in Kouveliotou 
et al. (2003b); an accurate (< 0".3) location of the source 
is given in Wachter et al. (2003). During this observation 
IGR 16358 - 4726 fell off axis on the (front illiuminated 
CCD) ACIS-S2. The source was very bright, however, 
the pile-up effect was negligible due to the significant 
broadening of the point spread function (PSF), enabling 
reliable spectral measurements (§3.2). We observed the 
source again (Obs2) on 2003 April 21 for 47 ks; we se- 
lected the continuous clocking (CC) mode to mitigate 
any pile-up and to exploit the mode's 2.85 ms time res- 
olution. 

3.1. Timing Analysis 

Strong pulsations with a period ~100 min are directly 
visible in the light curve (Figure 1; upper panel) dur- 
ing Obsl (Kouveliotou et al. 2003a). We measured the 
pulse period by first maximizing the signal-to-noise in 
the light curve by selecting photon energies in the range 
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Fig. 1 — Top panel: Chandra ACIS-S2 lightcurve (1.0-10.0 keV) 
of IGR 16358 - 4726 in an active state on 2003 March 24. The 
modulation at 5880(50) s is clearly seen. Bottom panel: Chandra 
ACIS-S3 lightcurve of the same source on 2003 April 21. Note the 
different duration of the x-axis in the two panels. 

2 — 10 keV and by correcting the arrival times to the So- 
lar System barycenter. Next, we constructed a template 
pulse profile from the full data set using a trial period 
determined from an epoch fold search. The data were 
split into 4 roughly equal time segments and folded on 
this trial period. We compared the relative phases of the 
four segments to the pulse template and finally we fit- 
ted them to a linear trend, which yielded the correction 
to our trial period. After a series of iterations, we find 
the best fit period for Obsl to be 5880(50) s. The pulse 
profile has a dominant sinusoidal component with sig- 
nificant substructure and varies as a function of energy 
(Figure 2; left panel). A portion of the substructure can 
be attributed to a quasi-periodic oscillation feature at 
~ 3 x 10~ 3 Hz. The pulse fraction is constant with en- 
ergy across the Chandra bandpass; the 2 — 10 keV rms 
value is 36.6(8)% (70(6)% peak-to-peak). 

During Obs2 the source was detected at a lower in- 
tensity, a factor of ~ 15.5 fainter in flux (2 — 10 keV). 
In addition, a visual inspection of the light curve shows 
substantially more flickering than present in Obsl (Fig- 
ure 1; lower panel). The Fourier power spectrum shows 
the 100 min pulsations at much weaker significance su- 
perposed upon an excess of red noise. There are no other 
coherent high frequency signals out to 175 Hz. Using 
the same technique as before, we measure a pulse period 
of 5860(30) s, consistent with the period derived from 
Obsl. The folded pulse profile during Obs2 (Figure 2; 
right panel) is significantly different than what it was 
four weeks earlier; the pulse shape is more complex with 
more power in the higher harmonics. Despite the changes 
in pulse shape, the pulse fraction has remained high and 
constant at 39(2)% rms (2 — 10 keV). We performed a 
cross-correlation of the high and medium energy (Figure 
2) pulse profiles of Obs2 and find a soft lag of ~ 110s; a 
lag was not detected in a preliminary analysis of Obsl. 
A refined analysis will be reported elsewhere (Patel et 
al. in preparation). 

3.2. Spectral analysis 

IGR 16358-4726 is detected significantly off-axis (9.7') 
during Obsl, necessitating a large extraction radius for 



ASCA/GIS 2-10 keV Flux a 

Date Exposure Time (ks) Absorbed Unabsorbed 

Feb 1999 69 11 24 

BeppoSAX MECS/LECS 2-10 keV Flux 
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Fig. 2.— Folded light curves of IGR 16358 - 4726 in an active 
state taken with the Chandra/ ACIS (in three energy bands) on 
2003 March 24 (Left paneJ) and 2003 April 21 (Right panel). 

the source spectrum to account for the significant PSF 
distortion of the source and for the scattering of the hard 
photons of the source spectrum. We have, thus, collected 
source events from a circular region of 1.5' radius. We 
extracted a background spectrum from a source free re- 
gion of 3.0' radius and on the same CCD. Finally, we 
grouped the data into energy channels that contained at 
least 25 events each. 

We report here results on our phase-averaged fits to the 
data. During Obsl the spectrum includes a highly sig- 
nificant Fe line feature and is hard and highly absorbed. 
We have initially fitted the data with a PL + a Gaussian 
line model; the best fit parameters are given in Table 
2, and the spectral fit is shown in Figure 3. The 2003 
Chandra flux is ~ 70 times higher than the ASCA flux 
of 1999. We have also fitted the data with an absorbed 
thermal blackbody (BB) model with a line (Table 2); the 
fit statistic is slightly improved. Finally, we have tested 
two-component (continuum) models of PL + BB and PL 
+ PL; we also find acceptable fits. However, the addi- 
tional (PL) component in each case is very soft and steep, 
mostly contributing to a very narrow spectral range of 
the continuum below 3 keV, where the signal is very low. 
We conclude that the single model fits (PL or BB) better 
describe the data. 

During Obs2, the spectrum is also well described by an 
absorbed PL. We find that the absorption column has de- 
creased significantly since Obsl, but is still in agreement 
with the best fit column density measured by ASCA, sug- 
gesting that at least part of the column is due to local 
absorption. Furthermore, we find that the Fe line feature 
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Fig. 3.— Chandra ACIS-S3 spectrum of IGR 16358 - 4726 in an 
active state on 24 March 2003 (top curve) and again on 21 April 
2003 (bottom curve). In both observations we find that a highly 
absorbed single power law model gives the best fit to the data. A 
Fe K a line feature clearly seen at 6.4 keV during the 24 March 
observation had faded below the detection capability of Chandra 
during the second observation. 

has faded beyond detection in Obs2. The upper limit we 
measure on the line equivalent width suggests that the 
flux in the line has faded more than the continuum flux; 
the latter is now a factor of 15.5 lower than Obsl. Fig- 
ure 3 also shows the PL spectral fit result of Obs2. The 
spectral parameters are listed in Table 2, together with 
the results of an absorbed BB model fit. 

Finally, we have fitted the two Chandra observations 
simultaneously using an absorbed PL + Fe line model. 
Here we linked the galactic absorption and Fe line com- 
ponent and added an additional intrinsic absorption com- 
ponent; the latter was free to vary for Obsl but assumed 
to be zero for Obs2. With 7 free, we derive galactic and 
intrinsic (Obsl) column densities and photon power law 
indices consistent with the individual fit values above 
(Table 2). In a similar fashion, we simultaneously fit 
both observations using an absorbed BB + Fe line model. 
With Nji t g a i and line parameters linked and UFbb free, 
we find galactic and intrinsic (Obsl) column densities 
and temperatures also consistent with the previous (not 
linked) values (Table 2). 

Further, we performed phase-resolved spectroscopy in 
five phase-bins with the data of Obsl using a PL + a 
Gaussian line model. We find that while the column 
density remains constant, the spectral index rises slightly 
during the ascending part of the pulse profile; the spec- 
trum is marginally softer during the pulse decay. The 
6.4 keV Fe line intensity is constant throughout the en- 
tire pulse. We will report detailed phase-resolved spec- 
troscopy results for both observations in a future work 
(Patel et al. in preparation). 

4. DISCUSSION 

Our Chandra observations of the Fe line centroid in the 
IGR 16358 — 4726 spectrum indicate that most likely the 
source is in a galactic X-ray binary system. If the 1.6 
hr period is orbital, then there are two types of X-ray 
binary systems which can have similar periods and X- 
ray light-curve properties: a Cyg X-3 like system (4.8 hr 
period; Milgrom 1976; Parsignault et al. 1977) and high- 
inclination accretion disk corona (ADC) systems such as 
4U 1822 - 37 (5.57 hr period; White et al. 1981). Like 
both these sources IGR 16358 — 4726 in outburst is show- 
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Table 2 

Chandra SPECTRAL FIT PARAMETERS" 



Model 


Tl 


r 2 


kT BB 


N H ,Gal iVff.snt E Hne EW 


F x (Abs) E. 


r (Unabs) 


^linc 


X 2 /DOF 


OBS1 ONLY 


PL+Gaussian 


0.5(1) 






3.3(1) ... 6.39(1) 0.14(6) 


6.9 


17.1 


0.30 


357.5/316 


BB+Gaussian 






3.1(1) 


2.90(8) ... 6.39(1) 0.13(3) 


6.7 


14.2 


0.27 


346.1/316 


PLi+PL^+Gaussian 6 


0.8(1) 


10(1) 




4.0(2) ... 6.38(1) 0.11(3) 


6.8 


66.2 


0.29 


337.6/314 


BB+PL 2 +Gaussian c 




8(1) 


2.8(2) 


3.5(2) ... 6.38(1) 0.11(3) 


6.6 


28.7 


0.26 


330.5/314 


OBS2 ONLY 


PL+Gaussian d 


0.8(2) 






2-0(1) 


0.6 


1.1 


< 0.02 


124.3/136 


BB+Gaussian d 






2.6(2) 


1-6(1) 


0.5 


0.9 


< 0.02 


123.5/316 


OBSl+OBS2 e 


PL+Gaussian d 


0.5(1), 0.7(2) 






2.0(1) 1.3(2) 6.39(1) 0.14(6) 


6.9,0.6 


17.1,1.1 


0.24,< 0.02 


478.8/448 


BB+Gaussian d 






3.1(2), 2.6(2) 


1.6(1) 1.3(1) 6.39(1) 0.13(3) 


6.7,0.5 


14.2,0.9 


0.21,< 0.02 


466.5/448 



Notes: Uncertainties are 68% confidence limits. All models include galactic absorption a kTBB, Ei ine and EW in keV, Nn in 10 23 
cm -2 , fluxes are in 10~ xl ergs s~ x cm~ 2 2 — 10 keV) b Double power law + Gaussian line model fits. PLi component constitutes 30% of 
the unabsorbed flux C PL2 corresponds to the additional power law model component in this model. BB component constitutes 56% of 
the unabsorbed flux d Line Flux: 90% flux upper limit for a line at 6.4 keV with similar intrinsic width as the one in Obsl e Data from 
Obsl and Obs2 are fit simultaneously. The pairs of values represent the spectral parameters from Obsl and Obs2 respectively. 



ing a smoothly varying X-ray light-curve; the depth of 
the modulation is not indicative of an eclipse and there 
is no phase dependent absorption. All three spectra are 
similar, well fit with a (hard) power law with a low en- 
ergy cut off at ~ 4 keV and with large, broad iron lines 
at ~ 6.4 keV. The source luminosities are also within 
the same range of ~ 10 36 erg/s (assuming a distance for 
IGR 16358 - 4726 of 10 kpc). 

It is, however, questionable if a He star with a strong 
stellar wind such as in Cyg X-3 could fit a system with 
an orbital period of 1.6 h. More likely, the system is an 
LMXB such as 4U 1822 — 37; the questions that immedi- 
ately arise are where did the extreme absorption detected 
in the spectrum of IGR 16358 — 4726 originate if not from 
the stellar wind of a massive companion, and is a low 
mass companion sufficient to produce the Fe line? One 
expects systems with such periods to be transient (King, 
Kolb & Burderi 1996) as the average mass transfer rate 
would probably be quite low. The 'outburst' accretion 
rate could nevertheless be high, thus, if the ADC inter- 
pretation is right, the true luminosity of the object is 
much higher than observed. It is just conceivable that 
this actually causes the intense absorption, if the accre- 
tion rate is highly supcr-Eddington for the duration of 
the outburst. This would produce a dense outflowing 
wind, as is probably seen in some of the Ultra Lumi- 
nous X-ray sources (Fabbiano et al. 2003), which could 
give the absorption. Similar super-Eddington outflows, 
with extremely high columns, are seen in quasars in- 
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ferred to be accreting at super-Eddington rates (Pounds 
et al. 2003; King & Pounds 2003). Sources like this or 
with even higher columns could easily have been missed 
by earlier satellites, so there could be a whole class of 
them. 

Finally, we would like to note that if IGR 16358 — 
4726 is a LMXB, its orbital period would be in the range 
of those of the ms pulsar subclass, only five of which 
have been found. No ms pulsations were detected so 
far and unless the source becomes active again, current 
instrumentation docs not afford ms pulsation searches. If 
on the other hand the 1.6 hr period is a neutron star spin 
period, the source is similar to 2S 0114+650, which has a 
spin period of ~ 2.7 hr (Finley, Belloni & Cassinelli 1992) 
and an Fe line at 6.4 keV (Yamauchi et al. 1990); by 
analogy to other slow accreting pulsars we would expect 
the magnetic field of the neutron star in IGR 16358—4726 
to be of order 10 12 Gauss. 
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